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boundary layer measured at 2m, 7m and 23m in Ris$, Denmark, 
were analyzed with particular emphasis placed on determining 
characteristics of the high frequency region of the spectra 
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one-dimensional temperature spectrum and an estimate of the 
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I. INTRODUCTION 



Several attempts over the past few years have been made 
to measure high frequency temperature variations in the 
atmospheric boundary layer. Sensors with a sufficiently 
high frequency response and spatial resolution have been 
devised to overcome the inherent problems previously 
encountered in trying to measure small scale, high frequency 
temperature fluctuations. In addition to the sensitive 
sensors, low-noise electronic equipment has become available 
to record and analyze small scale fluctuations. 

The object of this research was to utilize an analog 
analysis scheme to determine the shape of the high wave 
number (small scale) one- dimensional temperature spectrum 
and to evaluate the scalar constant, Kg ' , which appears in 
the expression for the -5/3 form of the temperature spectrum. 
In addition, simultaneous analyses of velocity fluctuations 
were conducted in order to obtain the necessary information 
ultimately required to evaluate Kg'. 

A knowledge of the small scale temperature distribution 
in a turbulent field has a number of applications. It is 
necessary when considering the heat budget which is of 
utmost importance in the fields of oceanography and 
meteorology. Sensible heat flux can be determined by 
the covariance of measured temperature fluctuations and 
vertical components of the simultaneous wind velocity 
fluctuations . 
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Temperature inhomogeneities in the atmosphere are 
one of the quantities associated with the scattering of 
acoustic waves and electromagnetic radiation. This 
scattering, which is a major concern in the fields of 
acoustics and optics, is often related to variations in 
the refractive index. 

Further applications of a knowledge of the small scale 
temperature distribution are found in the areas of physi- 
cal chemistry and thermodynamics. 
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II. THEORETICAL DEVELOPMENT 



A. VELOCITY SPECTRUM 

Small scale turbulence in a fluid is a feature of 
interest to the environmental scientist due to its effect 
on heat, momentum and energy transport. The one-dimensional 
velocity (kinetic energy) spectrum is given by 

u 2 = r $ (k) dk (1) 

0 

~2 

where u is the total variance of the turbulent velocity 
fluctuations and k is the radian wave number. The rate of 
energy dissipation in an isotropic, turbulent flow field is 
given by 



c = 15v f°° k 2 $ (k) dk (2) 

0 

2 

where the function k $(k), the dissipation spectrum, 

statistically describes the rate of decay of turbulent 

2 - 1 

kinetic energy to heat, v (cm .sec ) is the kinematic 
viscosity of the fluid. 

Modern turbulence research had as its impetus the 
concept of a continuous cascade of kinetic energy from 
large scale anistropic velocity fluctuations to smaller 
scales where the motion became increasingly isotropic and 
at which the energy is dissipated as heat. Small scale 
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motions which are characterized by small time scales are 
statistically independent of the relatively slow large 
scale turbulence derived from the mean flow. 

This concept was presented by the Russian mathematician, 
A. N. Kolmogorov in 1941. Kolmogorov proposed the idea of 
local isotropy, i.e., at some point in the cascade of 
kinetic energy from large scales to small scales, the 
turbulent fluid motion becomes isotropic. This idea was 
described by Kolmogorov in two hypotheses: 

(1) If the Reynolds Number (ratio of inertial to 
viscous forces) is sufficiently large, then there exists 
a range of small scales in which statistically, the turbu- 
lence is basically steady-state and is dependent only on the 
rate of dissipation of kinetic energy, e [ (cm/sec) ^ . sec ^)], 
the wave number, k(cm ^) and the kinematic viscosity, 
v(cm .sec ). Accordingly, dimensional analysis yields 

4>(k) = (ev 5 ) 1/4 F(k/k s ) (3) 

for the form of the one-dimensional energy spectrum and 

k 2 »(k) = k s W) 1/4 Ck/k s ) 2 F(k/k s ) (4) 

for the form of the energy dissipation spectrum. Here, 
k g = (e/v ) ' is the Kolmogorov wave number, l/k s is known 
as the Kolmogorov microscale. F(k/k s ) is a universal func- 
tion which is valid for all fields of turbulence provided 
the Reynolds Number is sufficiently large. 
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(2) If the Reynolds Number is very large (on the order 
of 10^ and greater) then a range of scales exists wherein 
the energy cascades independent of viscosity. In this 
range of wave numbers, the decay of the flow energy is 
characterized by a -5/3 dependence upon wave number. This 
range in which the energy is transferred inertially is 
known as the inertial subrange and characterizes the turbu- 
lence scales too large to be influenced by viscosity but 
small enough to be separate from the mean flow. In the 
inertial subrange, 



$(k) = K'e 2/3 k" 5/3 (5) 

where $(k) is the one-dimensional energy spectrum 
2 

[(cm/sec) .sec] and K' is a universal constant. Previous 
researchers have reported values of K' ranging from 0.48 
to 0.69 although the lower values are considered to be 
more accurate. 

B. TEMPERATURE SPECTRUM 

It is possible to apply the Kolmogorov ideas to 
turbulent fluctuations of passive scalars such as tempera- 
ture. Temperature is a passive scalar when it is trans- 
ported in a fluid but does not introduce buoyancy effects. 

If 0(C°) is the deviation of the temperature from its 
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mean, the one dimensional temperature fluctuation spectrum 
is defined as 



e 2 = f°° $ (k)dk (C 0 ) 2 . C6) 

0 0 

The rate of temperature dissipation in an isotropic flow 
field is 



e a = 6k r k 2 $ Q (k) dk CC 0 ) 2 .sec _1 . (7) 

0 0 0 

2 

k is the molecular diffusivity and k <J> Q (k) is the scalar 

dissipation spectrum which describes the distribution with 

2 

wave number of the rate of decay of the quantity 6 , the 
temperature variance. 

Again, according to Kolmogorov: 

(1) If the Reynolds number is sufficiently large there 
exists a range of small scale temperature fluctuations which 
statistically are dependent only on e, v, (the tempera- 
ture dissipation rate) and k. Once again dimensional 
analysis yields 

4> 0 (k) = e 0 e' 3/4 v 5/4 H(a,k/k s ) (8) 

for the one-dimensional temperature spectrum and 

k 2 * e (k) = k s 2 E e e‘ 3/4 v 5/4 (k/k s ) 2 II(o, k/k s ) (9) 
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for the temperature dissipation spectrum, a is the 
Prandtl number (a = v/k) , the ratio of kinematic viscosity 
to thermal diffusivity. In a gas a is approximately 1. 

For air at 20°C the value is near 0.7. 

(2) If the Reynolds Number is even larger, then the 
temperature spectrum is independent of v and k and is 
determined only by c and e 0 . Dimensional analysis 
[Corrsin, 1951] yields 

$ Q (k) = K 0 'e 0 e' 1/3 k" 5/3 . (10) 

K 0 ' is supposedly a universal constant although previous 
researchers have reported values ranging from 0.2 to 2.0. 

In view of this range of reported values of K 0 ' it 
was the intent of this research to determine the shape of 
one-dimensional temperature spectra at high wave numbers 
and to evaluate K 0 ' from the data analyzed. In this way, 
a larger data base could be established from which deter- 
minations for Kg’ could be resolved. 
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Ill . 



INSTRUMENTATION AND EXPERIMENTAL APPROACH 



A. BACKGROUND 

Sensors and electronic systems with the necessary 
sensitivity and response to measure high frequency, small 
scale atmospheric temperature fluctuations were developed 
in the early 1970's by National Electrolab Associates 
Limited of Vancouver, British Columbia, in consultation 
with the Institute of Oceanography, University of British 
Columbia. Previous electronics allowed the temperature 
spectrum to be measured to just beyond the inertial sub- 
range but not to sufficiently high wave numbers to 
unequivocally evaluate the area under the scalar dissipa- 
tion spectrum. 

B. TEMPERATURE MEASURING SYSTEM 

1 . Sensor 

The sensor of temperature fluctuations utilized 
was a resistance thermometer consisting of a platinum wire 
0.25 pm (0.00001 in.) in diameter and 0.30 mm in length. 

The wire has a time constant of less than 10 ys in a 4 m/sec 
air flow. With a detection current of 50 yamps, its 
sensitivity is 0.063 mv/°C. The sensor is made from 
Wollaston wire consisting of a silver jacket about a 
platinum core. The outside diameter of the jacket is 
45 ym (0.0018 in.). The jacket and core are prestressed 
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into a V- shape and then soldered onto the end of a standard 
hot wire anemometer probe (Figure 1) . The sensor has a 
frequency response which is essentially flat to 5KHz. 

2 . Bridge 

Utilizing a thin platinum wire to measure temperature 
fluctuations takes advantage of the fact that a slightly 
heated wire is basically a simple resistance thermometer. 

In order to sense a change in resistance in the wire due 
to variations in the ambient temperature, a small heating 
current is passed through the wire and the resulting voltage 
drop can be measured. If the heating is kept low, the 
resistance changes due to cooling by wind are negligible 
compared to resistance changes caused by temperature fluc- 
tuations. A bridge circuit is utilized to ensure optimum 
current flow through the platinum core. 

The bridge circuit is a modified Wheatstone bridge 
consisting of two 2.21Kft resistors, sensor, probe and cable, 
balance adjustment resistors and capacitors. The fixed 
and variable capacitors (Figure 2) paralleling the balance 
resistance are necessary to match the probe and cable 
capacitance. The bridge amplifier is a differential 
amplifier with a fixed gain of 60dB and a frequency 
response flat to 150KHz. 
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Figure 1. Wollaston Wire Mounted on Probe 
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Figure 2. Bridge Circuit 



C. EXPERIMENTAL PROCEDURE 



1 . Location 

The data were collected at the meteorological 
field site operated by the Meteorological Group of the 
Danish Atomic Energy Commission, Research Establishment, 
Ris?5. The site features a 130M meteorological tower 
instrumented at seven levels. In addition, miscellaneous 
instrumentation can be mounted at arbitrary levels provided 
only that the electronics is sufficient to drive the cables 
from the instrumentation to the junction boxes at the 
seven standard levels. The tower stands on the Ris0 
Peninsula 0.5 Km from Roskilde Fjord. Signals from tower 
instrumentation are carried by cable to a meteorological 
station 50m downwind of the tower where both analog and 
digital recording equipment are housed. During the course 
of the experiments winds blew essentially from the west 
off the fjord. Data were not collected unless the wind 
was from off the water. 

The Ris^ Peninsula is 6Km north of the ancient 
Danish town of Roskilde which in turn is 30Km approximately 
southeast of Copenhagen on the island of Sjaelland. 

2 . Equipment 

a. Meteorological Tower Arrangement 

Mean wind speed and direction, temperature 
and humidity were measured at 7, 23, 39, 56, 72 and 96m 
(Figure 3) . Sonic anemometer measurements were made at 
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1. Gate-house and fire 
station 

2 . Administration 

3. Accelerator Department 

4. Agricultural Department 

5. Workshop 

6. Phxsics Department 

7. Electronics Department 

8. Reactor Phssics 
Department 

9. Chemistry Department 

10. Lecture Hall 

11. Canteen 



12. Library 

13. Health Physics 
Department 

14. Engineering Department 
and drawing offices 

15. Helium plant 

16. Service and maintenance 

17. Research Reactor DR 1 

18. Waste treatment plant 

19. Metallurgy Department 

20. Hot Cells 

21. Research Reactor DR 3 

22. Meteoro lo gy station 



23. Isotope Laboratory 

24. Research Reactor DR 2 

25. Reactor Engineering 
Department 

26. 'Tandem Accelerator 
Department of the 
Niels Rohr Institute 

27. Kindergarten 

28. Guest-house 

29. Staff dwellings 

30. “Svaleholm" farmhouse 

31. “D> sku.-rgard" farmhouse 



Figure 3. Ris0 Research Establishment 
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2, 23 and 72m. Fine structure turbulence measurements 
were made at different heights during the course of the 
experiment. In general. Naval Postgraduate School (NPS) 
equipment was deployed at 2m, 7m or 23m. Oregon State 
University (OSU) equipment was always at 23m. University 
of California, San Diego (UCSD) equipment was placed very 
near the top of the mast at 123m. Danish Atomic Energy 
Commission (AEK) equipment was located at 2m. Each group 
deployed sensors for measuring turbulent temperature and 
velocity fluctuations. The main purpose of the experiment 
was to make simultaneous measurements of atmospheric 
temperature fluctuations by different systems, at the same 
level and at different levels. The temperature sensor in 
each case was a platinum wire. The diameter of the NPS 
wire was 0.25 ym, of OSU 0.625 ym and of UCSD both 0.25 
and 0.625 ym. By making these simultaneous measurements 
with similar sensors but different electronics, it was 
hoped (1) to find in what way they differed and (2) more 
importantly, to resolve some of the discrepancies surround- 
ing the value of the scalar constant Kg'. 

The main comparisons made in this thesis are 
between the NPS and OSU systems. 

b. NPS Platinum Resistance Thermometer 

At the 2m height the sensor was mounted in an 
aluminum holder which could hold as many as four such 
sensors. For the data discussed in this thesis the 
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temperature sensor was placed 5mm below a Disa hot wire 
anemometer. The other two positions were left vacant. 

At 7m the same aluminum holder was used but a 
Thermosystems probe was used for the hot wire. This probe 
was vertically aligned and then rotated to be within 5mm 
to the side of the temperature sensor. 

A somewhat similar arrangement was made at 
23m except it was more difficult to get the temperature 
probe close to the hot wire probe. There was a separation 
of approximately 1cm. 

c. Hot Wire Anemometer 

A variety of hot wire anemometers was used. 

They were Disa (AEK, UCSD) , Thermosystems (NPS, OSU) , 
and a new wind- vane probe system (AEK) designed by Larsen 
and Busch (1974) . 

Disa probes were mounted in the aluminum holder 
previously described whereas the other systems had indi- 
vidual clamps. 

3 . Recording Systems 

NPS, OSU and UCSD used analog tape recorders. 

AEK recorded digitally on magnetic tape. The NPS data 
analyzed for this thesis were recorded on an Ampex FR 1300 
tape recorder at 7 1/2 ips using FM electronics. At this 
speed the 3dB point is at 2.5KHz. Each NPS temperature 
signal was recorded both direct and differentiated (with 
respect to time) . The object of the differentiation is 
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to increase the level of the signal at high frequencies 
above the noise level. The differentiator circuit had a 
gain of 6dB/octave over the differentiating band of 
frequencies. This band could be varied from DC to 500Hz, 
1000Hz, 1500Hz, 2000Hz or 3000Hz. The frequency at which 
unity gain occurred could also be varied. Beyond the band 
of differentiation, the circuit rolled off at 6dB/octave. 
No other filtering of the signals was done. NPS hot wire 
signals were recorded directly without filtering or 
differentiation . 

The input level of the tape recorder was set to 
handle 1 volt rms signals. 
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IV. ANALYSIS PROCEDURES 



A. SELECTION OF DATA 

Generally when concerned with the analysis of high 
frequency turbulent temperature and velocity signals, 
it is not necessary to examine long sections of data. 

It was determined, however, in the preliminary analysis 
of this data that in order to obtain characteristic 
spectra, it would be necessary to extend the record length 
to 10 minutes. This was due to varying levels of the 
recorded signals, necessitating the utilization of a 
broadband amplifier to maintain the signal level at an 
almost constant level. 

The primary considerations for selecting a section of 
data for analysis were first that the temperature and 
velocity fluctuations be considered typical of the record 
and secondly, if possible, correspond to sections that were 
being analyzed at Oregon State University for future com- 
parison of results (Figure 4) . By trying to adhere to 
these considerations, the following data sections were 
analyzed: a 10 minute length of record from Run 5 commenc- 

ing at 1537 local Ris0 time on 30 August 1974, and a 10 
minute length of record from Run 6 commencing at 1657 on 
30 August 1974 (Table I) . 
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Figure 4. Example of Data Selected for Analysis 
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B. ANALOG EQUIPMENT AND PROCEDURES 



The analog analysis scheme consisted of utilizing the 
following state-of-the-art equipment to obtain spectra of 
the selected sections (Figure 5) : 

A Honeywell Model Ninety-Six Magnetic Tape Recorder/ 
Reproducer System was used to playback the recorded 
signals. This device is capable of recording/reproducing 
any combination of multi-channel direct (analog) and FM 
data at nine selectable servo-controlled tape speeds 
from 15/16 through 240 inches per second (ips) . Its 
standard magnetic assemblies are IRIG compatible. It is 
characterized by a tape speed accuracy of 0.1% and a 
bandwidth of 0-2.5KHz (within ldB) at a signal-to-noise 
ratio of 50dB when operated at 7-1/2 ips. 

The playback signals from the reproducing system were 
fed to a Preston 8300 XWB Amplifier with a bandwidth of 
lOOKHz and selectable gains of 1, 2, 5, 10, 20 and 50. 

The output from the amplifier was fed into the input 
of a Federal Scientific UA-500 Ubiquitous Spectrum Analyzer. 
The selected analysis range was 0-2KHz at a sample rate 
of 6000hz. The analyzer has a 3dB bandwidth of 6hz and 
a frequency accuracy of ±0.2% of the full analysis range. 

The input signal to the analyzer was monitored by a 
Hewlett-Packard 3400A RMS Responding Voltmeter and a 
Tecktronix Type 502A Dual-Beam Oscilloscope. 
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Figure 5. Analog Analysis Scheme 



The spectral output was visually monitored on a 
HP 1300A X-Y Display and plotted on a Mosely 7035B X-Y 
Solid State Recorder produced by Hewlett-Packard. 

C. DATA INTERPRETATION 

Theoretical work almost always expresses turbulence 
results in terms of wave number while actual turbulence 
measurements are spectrally computed as functions of 
frequency. These two approaches are related by Taylor's 
hypothesis 



k = 



2irf 

U 



(ID 



where k is radian wave number (cm j , f is frequency (Hz) 
and U is mean wind speed (cm/sec). Equation 11 implies 
that turbulence with scales near 1/k are "frozen" during 
the period these scales are swept past the sensor. The 
maximum frequency analyzed was 2KHz at a mean wind speed 
of 5 m/sec which indicated that a maximum wave number of 
25cm 1 would have to be resolved. This requires that the 
sensor does not exceed 0.4mm in length and have a fre- 
quency response essentially flat to 2KHz. Both of these 
requirements were met. 

Since Eq. 11 is linear, the shape of a wave number 
spectrum is the same as that of a frequency spectrum. 
Frequency spectra and wave number spectra are related by 
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= k<J> Ck) . 



( 12 ) 



The inertial subrange expression (eq. 5) in terms of f 
becomes 

<j>(f) = (g 2 -) 2/ K'e 2/3 f' 5/3 . (13) 

The rate of kinetic energy dissipation expression (eq. 2) 
becomes 

2 

e = lSv (jy2-) I™ f 2 <j>(f)df • (14) 

u 0 

The resulting expressions for the . temperature spectrum in 
terms of f are 

♦ e tf) - #f 2/s w' 1/s f ' 5/3 (15) 

and 

2 

e 0 = 6 k C^L) o /“f 2 <J) 9 (f)df . (16) 

Equations 13, 14, 15 and 16 were the expressions used 
to evaluate e, e Q and K A ' . 
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V. RESULTS 



A. VELOCITY SPECTRUM RESULTS 
1 . Calculation of e 

a. Indirect Method 

An indirect estimate of e was made from the 
friction velocity and the logarithmic wind profile 
(Figure 6) . Thus , 

Z + Z 

U* = KU/ J£N — - - (17) 

L 0 

where K = 0.4 is the von Karman constant and Zq the rough- 
ness length. The relation 

e = U* 3 /KZ (18) 

was then used to evaluate e . 

b. Inertial Range Method 

Once the spectrum 4> ( f ) was determined, e was 
estimated by choosing values of <j>(f) from the best straight 
line in the frequency region in which a -5/3 form was 
observed and then substituting into Eq. 13. A value of 
0.5 was assumed for the constant K' based on results 
presented by Boston (1970) and Williams (1974) . 

c. Direct Method 

The kinetic energy dissipation rate was 
estimated directly by summing of the frequency range 
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U (m/sec) 



Figure 6. Wind Profiles Runs 5 and 6 
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from 10Hz to 2KHz the scales contributing to viscous dissi- 
pation according to Eq. 14. 

Table II summarizes the results obtained 
from these three methods. All values of e determined by 
the three methods were felt to be representative of the 
actual value. Therefore an average value for each run 
at a given height was utilized in the further analysis 
of the temperature spectrum. 

2 . Normalization of Spectral Results 

a. Velocity Spectra 

The velocity spectra were normalized according 
to Eq. 3 so that a comparison could be made between dif- 
ferent runs and previous experiments. Figures 7 and 8 
are normalized velocity spectra from Runs 5 and 6 at the 
2m height. 

These figures show that the spectra nearly 
overlap with a characteristic break from the -5/3 region 
near log k/k g = -1. 

b. Energy Dissipation Spectra 

The energy dissipation spectra were normalized 
according to Eq. 4. A linear plot was used to display 
results (Figure 9) . Maximum dissipation occurred near 
the value of k/k s = 0.1 which is in agreement with previous 
research [Boston, 1970]. 
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COMPARISON OF e VALUES 
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Figure 7. Normalized NPS Velocity Spectrum Run 5 at 2 m 
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Figure 8. Normalized NPS Velocity Spectrum Run 6 at 2m 
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Figure 9. Normalized Energy Dissipation Spectrum 
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B. TEMPERATURE SPECTRUM RESULTS 



1. Calculations of e 0 

e Q values were estimated from Eq. 16 and the 

O 

results are included in Table III. 

2. Calculations of K Q ' 

u 

The values of K Q ' were determined from Eq. 15. 
Values of $g(f) were chosen from the best straight line 
in the frequency region in which a -5/3 form was observed. 
The value of e used was the average value for a given run 
at a particular height. Since e enters into Eq.. 15 to the 
-1/3 power, errors in estimating e do not seriously affect 

V* 

The results of the Kg' calculations are shown in 
Table III. An average value of 0.88 was obtained which is 
in quite good agreement with 0.81 by Boston (1970). It 
was noted that although no inherent differences in tempera- 
ture spectral shape were observed between OSU and NPS 
sensor results, the values of K ' determined at the 23m 
height from the OSU sensors were very near the values 
reported by Williams (1974) at OSU. 

3 . Normalization of Spectral Results 
a. Temperature Spectra 

The temperature spectra were normalized 
according to Eq. 8. The temperature spectra from each 
run at each height were normalized and plotted (Figures 
10-15). A composite spectrum was obtained by superimposing 
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Figure 10. Normalized NPS Temperature Spectrum Run 5 at 2m 
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Figure 11. Normalized NFS Temperature Spectrum Run 6 at 2m 
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Figure 12. Normalized NPS Temperature Spectrum Run 5 at 7m 
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Figure 13. Normalized NPS Temperature Spectrum 

Run .6 at 7m 
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Figure 14. Normalized OSU Temperature Spectrum 

Run 5 at 23m 



45 



00 k 



:> 

to 

tn 



CJ 



® 

•©- 



o 



bO 

O 

-J 



3 



2 



1 



0 



-1 



-2 



-3 





0 




O 




O 




O 




G 




G 


- 


G 




G 




O 


- 


G 


- 


G 

1 




O 

l i 1 



-3 -2 -1 0 



L °SlO k/k s 

Figure 15. Normalized OSU Temperature Spectrum 

Run 6 at 23m 
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ONE-DIMENSIONAL SCALAR CONSTANT K 
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23m 7.16 (6) 1 . 40 ( - 5 ) 0.00048 0.948 



all the spectra on a single plot (Figure 16) . The break 
from the -5/3 region occurs just beyond log k/k s = -1. 

These spectra show the shape of the one-dimensional 
temperature spectra in air beyond the -5/3 region which 
was one of the primary objectives of this thesis. The 
slope of the spectrum beyond the -5/3 break tended to be 
-32/5 which agrees well with Heisenberg (1948) who demon- 
strated that a slope of -7 should exist for very large 
values of k. With the exception of the temperature 
spectra at 2m, there is very little scatter in the data 
at high wave numbers . 

Figure 17 is a composite plot of characteristic 
temperature and velocity spectra. It is noted that there 
is very little difference between the two spectra but the 
temperature spectrum extends to slightly higher wave 
numbers than the velocity spectrum. A -32/5 slope is 
observed for both temperature and velocity. 

b. Temperature Dissipation Spectra 

The temperature dissipation spectra were 
normalized according to Eq. 9 and again plotted linearly. 
Figure 18 shows a typical normalized temperature dissipa- 
tion spectrum. There is very little difference between 
the temperature and velocity dissipation spectra. The 
normalized shapes are almost identical with the temperature 
spectrum dropping off slightly faster than the velocity 
spectrum . 
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Figure 17. Comparison of Normalized Temperature and 

Velocity Spectra 
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Figure 18. Normalized Temperature Dissipation Spectrum 
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VI. SUMMARY AND CONCLUSIONS 



The two primary objectives of this research were to 
determine the shape of the high wave number one -dimensional 
temperature spectrum and the evaluation of the Kolmogorov 
scalar constant Kg'. 

As seen in the results section, this investigator has 
determined that the high wave number one-dimensional 
temperature spectrum falls off with a slope of -32/5 beyond 
the -5/3 break which occurs slightly beyond log k/k s = -1. 
Also the determination of the Kolmogorov scalar constant. 
Kg', to be 0.88 is considered to be significant not only 
in absolute value but also that it increases the data 
base of high frequency temperature fluctuation analysis. 
Further, the criticism that previous evaluations were based 
on records that were too short is removed. It was also 
observed (Table III) that K Q ' increased in value both with 
height and Reynolds Number. 

As an additional result of this research, it is con- 
cluded that the spectra of temperature and velocity at 
high wave numbers are very similar with the temperature 
spectra extending to slightly higher wave numbers than 
the velocity spectra but that they both fall off at the 
same rate. 
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